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Abstract 

We propose a new model with the Majorana neutrino masses generated at two-loop level, in 
which the lepton number violation (LNV) processes, such as neutrinoless double beta decays, 
are mainly induced by the dimension-7 LNV effective operator O 7 = /^ 7 ^Li(D^<l>)<l> 4 >. Note 
that it is necessary to impose an Z 2 symmetry in order that O 7 dominates over the conventional 
dimension-5 Weinberg operator, which naturally results in a stable Z 2 -odd neutral particle to 
be the cold dark matter candidate. More interestingly, due to the non-trivial dependence of the 
charged lepton masses, the model predicts the neutrino mass matrix to be in the form of the normal 
hierarchy. We also focus on a specific parameter region of great phenomenological interests, such as 
electroweak precision tests, dark matter direct searches along with its relic abundance, and lepton 
flavor violation processes. 
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I. INTRODUCTION 


The presence of the tiny neutrino masses and mixings between different neutrino fla¬ 
vors have been established by many neutrino oscillation experiments l|-l7|, while more and 
more evidences are accumulated for the existence of dark matter (DM) over the last sev¬ 
eral decades, with the most precise measurement of its relic abundance by PLANCK 8|, l9|. 
Both phenomena cannot be explained within the Standard Model (SM), thus providing us 
with two windows towards new physics beyond it. An interesting idea is to connect neu¬ 


trinos anc 
Refs. 10 


DM in a unified framework as many existing attempts in the literature (see e.g. 


13|L We would like to push this connection further in the present paper. 


In order to understand the mass hierarchy problem in the neutrino sector, there are 
many models in the literature to naturally generate small Majorana neutrino masses such 
as the traditional Seesaw 


14 


26j and radiative mass generation mechanisms 


10 


13 


27 


21 . 


Most of them can be summarized as a specihc realization of the conventional dimension-5 
Weinberg operator. However, the generation of Majorana neutrino masses only requires the 
lepton number violation (LNV) by two units, and there exist many other equally legitimate 


LNV effective operators 


30 


, which are composed of the SM fields but with higher scaling 


dimensions. From the effective field theory perspective, it is generically believed that these 
high-dimensional effective operators are subdominated by the Weinberg operator due to 
the suppression from the corresponding high powers of the large cutoff. In order for these 
operators to show up as the leading contributions, one usually needs to impose an additional 
symmetry on the model to break the usual scaling arguments. Furthermore, if this symmetry 
is kept unbroken, then the lightest symmetry-protected neutral particle would provide a 
perfect DM candidate. In this way, the symmetry connects Majorana neutrino masses and 
DM physics by the high-dimensional effective operators. Such a connection has been already 


exemplihed by some recent three-loop neutrino mass models ^ 39-41], which realize the 
dimension-9 effective operator Og = with the DM embedded in the loop. 

In this study, we focus on a specific dimension-7 LNV operator Oj = 


32 


34 


37| and construct a UV complete model with an unbroken Z 2 symmetry to accomplish 


the above general arguments. In this model, Majorana neutrino masses arise radiatively at 
two-loop level, and neutrinoless double beta decays are dominated by a new “long- 
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FIG. 1: Typical diagrams for (a) induced neutrino masses by the effective operator O 7 , (b) the 
one-loop realization of O 7 , and (c) the two-loop neutrino mass generation. 

range” contribution as the results of the existence of Oj, while DM can also be embedded 
naturally as the lightest Z 2 -odd neutral state. 

This paper is organized as follows. In Sec. 11, we hrst describe the particle content and 
write down the relevant part of the Lagrangian for the model. We then calculate the two- 
loop neutrino masses and new contributions to the decay rate in the model, with 

the emphasis on their relations to the high-dimensional effective operator O^. In Sec. Ill, 
the constraints on the model are addressed from the electroweak (EW) precision tests, dark 
matter searches, and lepton flavor violating (LEV) processes. Finally, we give the conclusions 
in Sec. IV. 


II. GENERATION OF NEUTRINO MASSES AND On/3/3 DECAYS 


A. The Model 


Fig. [It shows the neutrino mass generation induced by the one-loop diagram with O 7 . 
In order to produce C >7 at one-loop level, we introduce two scalars: s : ( 1 , 2 ) and y : 
(2,1/2), and three vector-like fermions : (2,1/2) with i = 1, 2 and 3 to the SM under 

SU(2)i X U(l)y. K Z 2 symmetry is also imposed, in which only the new particles carry odd 


^ The definitions of “short” and “long” range contributions to On/3/3 follow Refs. 
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charges. The relevant new parts of the Lagrangian are given by 


— AL = s) + x) + + ^ii{DLiY{ia 2 )Lus + Cii{DLi)lmX + ^DiDLiDiii 


+KS*x\i(T 2 )^ + y + H.c. 

A3(<h^<h)(x^x) + A4(<htx)(x^<h) + A$.(<h^<h)(s*s) + A;,,(xtx)(s*s) 
+-^x(x^x)^ + As(s*s)^ , 


1/4 = 


( 1 ) 

( 2 ) 


where <h = (‘h^, is the SM Higgs donblet and a 2 is the Panli matrix for the SU(2)j;, 
gange gronp. After the EW spontaneons symmetry breaking, $ acqnires a vacnnm ex¬ 
pectation valne v = \/2{<h‘^) ~ 246 GeV, while > 0 is necessary for preserving the Z 2 
symmetry. Notice that the lepton nnmber is explicitly broken only when k, A 5 , and at least 
one of are non-zero simnltaneonsly. For convenience, we dehne + |A$sn^ 

and + The trilinear conpling constant k, in Eq. ([T]) makes mix with y;=*=, 

which can be formnlated as 




X 


± 



Sf \ \/2 kv 

S± ) ’ ^ M2 - Ml ’ 


( 3 ) 


where = sinx, = cosx, = tanx, and and represent the charged mass 
eigenstates with their masses, given by 


Mi = 


-Mi 


C2e 


C 29 


-Ml. 


Ml = - ^Ml , 

029 ^ 029 


( 4 ) 


respectively. On the other hand, A 5 contribntes to the mass splitting between and A° in 
X, shown as 


Mjj — M'^ -\- 2 (-^4 + As)^^ , M^ — M^ -|- -(A4 — As)^^ . 


( 5 ) 


As for the new fermions, we have the tree-level relation M^± = Mj^o for each Di donblet. 
The mass splittings between the charge and nentral components of the inert fermion donblets 
can only be indnced by loop corrections with valnes aronnd a few hnndred MeV j^. 

In this paper, we will characterize the model by nsing the physical qnantities: 


Mh. Ma. Ms^.Ms 2. MJ:)^, se, Xl , in. and On 


( 6 ) 
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where \l = |(A 3 + A 4 + A 5 ), and the other independent coupling constants from quarter 
terms: 


A$S5 Aj^S; As, 

which are less relevant in our discussion. 


(7) 


B. Two-Loop Majorana Neutrino Masses 


As seen in Fig. [T]d, the effective operator 0^ can be induced by the one-loop diagram, 
whereas the Weinberg operator cannot^. Consequently, Majorana neutrino masses appear 
through the two-loop diagram in Fig. The resulting neutrino mass matrix My dehned 


in the Lagrangian —+ H.c. can be calculated as 

1 GpS2e 


{My)iii = 


(167r2)2 ^ 


{Mjj — M^) ^2 Cil ill' + Cil' ^il) [lil — 7 * 2 ] , 

i=l,2,3 


where mi (/ = e, p, r) are charged lepton masses, and lij are dehned by 


lij — 


i>l rl—j/2 r^—^3 pl—x2—x3 

/ dy2 / dyi / dx^ / dx2 / dxi 

'0 JQ Jo Jo Jo 


2(1-3a:) 67 / 1 ( 2 -x) 

a;(l — x) (1 — xy 


iogH) + 

a:(l — X) mfj 


mfj = yi[xiMfj + X2M^ + xsM^] + y2x{l - x)M2 
+ - yi - y2)x{l - x)Ml^ , 

X = Xi+ X2 ■ 


( 8 ) 


(9) 


( 10 ) 


Subsequently, one can diagonalize My by 


diag(mi,m2,m3) = V'^ MyV, 


( 11 ) 


where mi^ 2,3 are three neutrino mass eigenvalues, which can have the normal ordering, mi < 
m 2 <C m 3 , or inverted ordering, m 3 -C mi < m 2 , and V is the Pontecorvo-Maki-Nakagawa- 
Sakata mixing matrix j^, ^|. Without loss of generality, V can be written as the standard 


^ There is a similar realization of O 7 in Ref. H, in which a triplet replaced the singlet s of our model. 
A fundamental distinction of their paper from the present one is that (P 7 does not give the dominant 
contribution to Majorana neutrino masses in Ref. 371. 

^ Similar topology with one exchange in a two-loop neutrino mass model can also be found in Refs. 

46j| . in which a different high-dimensional effective operator is realized without DM. 


45 
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parametrization by appropriate rephasing in L^’s and Irs, given by j^| 



C12C13 S12C13 5136"*'^ 


0 

0 

V = 

— ■S12C23 — Ci2S23'Sl3e*^ C12C23 — Si2S23Sl3e*'^ S23C13 


0 e*"2i/2 0 


i "S12S23 — Ci2C23Si3e*'^ —C12S23 — Si2C23Si3e*^ C23C13 , 


0 0 


.( 12 ) 


where the mixing angles with Sij = sinOij and Cij = cos 6ij are dehned within [ 0 , vr/ 2 ], and 
Dirac phase 6 and Majorana phases a 2 i and a^i are dehned within [0,27r]. 

From Eq. (|H]), we can get two important features for this mass generation mechanism. 
Firstly, the overall size of My is proportional to the mass difference of the neutral scalars, 
Mjj — M\, and the combined factor of the charged states, S 2 e{Iii — Ii 2 ), in which the former 
is generated by As and the latter corresponds to the size of k. Turning off one of them will 
make all neutrinos massless. Secondly, the neutrino masses are positively correlated to the 
coupling matrix elements and (u, as well as the sizes of m/. As the existence of the 
charged lepton mass hierarchy, nie -C < nir, if both matrices of {^u) and {(u) are in 
uniform textures, the magnitude of (M^)ee should be much smaller than those of other My 
entries. We make a great advantage of this general expectation by taking the following limit 

{My),, ~ 0 , (13) 


which is shown in Refs. 


49, 


50| to rule out the inverted ordering of neutrino mass spectrum at 


more than 2cr conhdent level. Thus, the normal ordering is predicted for the present model. 
Note that in the limit of Eq. flT5]l . Ref. ^ even shows that the lightest neutrino mass mi 
can only be located within the range O.OOleV < mi < O.OleV. Moreover, the smallness of 
{My),, is also required by the 0i>f3[3 decay, which will be clear in the next subsection. 

If we further focus on the CP conserving case, i.e., 6, a 2 i, 031 = 0 or vr, then mi will 
be constrained in the two narrow regimes, around 0.002 and 0.007 eV with {5, 021 , 031 } = 
{0, TT, 0} (Texture A) or {tt, tt, 0} (Texture B) and {0, tt, tt} (Texture C) or {tt, tt, tt} (Texture 
D), respectively. Taking the central values of ^ 12 , 023, 0i3, Am^i, and Am 32 from the global 


htting for the neutrino oscillation data 


49| . the corresponding mass matrices for Textures 


A, B, C, and D (Ta, Tb, Tc and Td) are given by 


^ 0 0.12 0.92 ^ 


Ta : My — 


0.12 1.9 2.7 
0.92 2.7 2.4 


(10-^) eV, 


(14) 


/ 
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^ 0 - 0.90 - 0.24 ^ 


Tr : M. = 


Tc : M, = 


/ 


Td : My 


- 0.90 

1.7 

2.7 

- 0.24 

2.7 

2.6 

0 

- 1.1 - 

0.055 

- 1.1 

- 2.3 • 

- 2.1 

0.055 

- 2.1 • 

- 3.1 

0 

- 0.086 

1.1 

■ 0.086 

- 2.6 

- 2.2 

1.1 

- 2.2 

- 2.9 


(10-^) eV, 




(10-^) eV, 




(10-^) eV, 


( 15 ) 


( 16 ) 


( 17 ) 


respectively. 

We now search for possible coupling matrix forms to realize the above four CP conserving 
neutrino textures. For simplicity, hereafter we will take M/ji = M £,2 = = Md, and also 

set ^ proportional to the identity matrix with the diagonal matrix element to be ^d- Taking 
a symmetric form of (, the mass matrix element should be proportional to + t^v)- 

We remark that by an appropriate phase absorption to the fermion helds, one can always 
have a positive ^d without loss of generality. Comparing with Eqs. ffTT)) - ffT7)l . the forms of 
Qii! in the four neutrino matrix textures can be obtained as 



X 

0.12 0 . 052 ^ 



X 


- 0.84 

- 0 . 013 ^ 


Ta : C oc 

0.12 

0.89 0.14 

, Tb : 

C oc 

- 0.84 

0.82 

0.14 



^ 0.052 

0.14 0 . 068 ^ 



^ - 0.013 

0.14 

0.072 ^ 



X 

- 1 . - 0.0031 ^ 



/ 

X 

- 0.081 

0.062 ^ 

Tc : C « 

- 1 . 

- 1.1 - 

- 0.11 

, Td 

: C oc 


- 0.081 

- 1.2 

- 0.12 


^- 0.0031 - 0.11 - 

0.088 ^ 



V 

0.062 

- 0.12 

- 0.081 ^ 


( 18 ) 


where the cross means that the value of C,ee is still arbitrary at this stage, which will be 
constrained by Of(3(3 decays. The overall scale of (hi can be determined by Eq. ([S]) when 
all new particle masses, sg and ^d are known. We will also leave the discussion about the 
correlation between ^ and ( from the LEV constraints to Sec. III. 
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Ml 


Ml 





(b) 


(c) 


FIG. 2: Diagrams for Oz^/3/3 decays from (a) the traditional processes with the neutrino mass inser¬ 
tion on the propagator, (b) the contribution involving O 7 , and (c) the one-loop construction which 
realizes Oy. For (b) and (c), the corresponding upside-down diagrams also need be considered. 

C. Neutrinoless Double Beta Decay 


In the previous section, we have built a two-loop neutrino mass model in which the LNV 
operator Oy provides the leading contribution. The next step is to study the LNV effect in 
this model induced by this high-dimensional operator. The most sensitive smoking gun for 
the LNV is the 0z//9/3 decay, which will place the strong constraint on ^d(ee and 

For the Majorana neutrino masses, there always exists the traditional long-range decay 
process by the exchange of neutrinos with a pair of left-handed electrons emitted as shown 
in Fig. [2^. Note that there is a chirality flipping on the internal neutrino propagator, which 
leads to the proportionality of the amplitude to the neutrino mass matrix element {M^)ee, 
in the sense that the detection for Oi//3/3 processes could help to determine or constrain 
\{M^)ee\. However, it is known that in some types of neutrino models. Fig. [2 |l does not give 
the main part of this process, and one should have a prior consideration on the effects of other 
new diagrams. For example, a class of neutrino models 


51 


53| that can be characterized 
by the dimension-9 operator Og = r^/i^[(Zi)^<F)'^(fcr 2 )<F][(Tl^<F)^(fcr 2 )<F] is well studied in 


Refs. 


51 


52 


54 


5fil|. and the new contribution is much larger than that from Fig. [2^ by 
orders of magnitude of 10®. For our model or those with Oy as the main LNV source, 
the 0z//5/3 decays are dominated by the diagram in Fig. [2 ]d, which is not suppressed by the 
nearly-vanishing (M,^)ee. We can write down a general formula for the half lifetime Ty 2 ^ of 

























the Oi//3/9 decay with the contributions from Figs. [2^ and[2]D, given by |h 




V = 


me 

1 S2eiC 


me 


IGtt^ 




(19) 

( 20 ) 


with 


»l-a:3 


'‘I-X2-X3 


= 




dx-} 


dxi 


XiMfj + X 2 M\ + X3MI. + (1 - Xi - 0)2 - x^)M\ 


D ’ 


where Cmm, Cm 
dehned in Ref. 


^and Cr^r] include the phase space integrations and nuclear matrix elements 


57l |. and r] is the coupling of the interaction {AGp/\/2){uL'y^dL){lR'^^vj}-i 
OOP generation for O-j (in Fig. [2J::). By using the numerical 
, we hnd that the contribution proportional to Cr^ri is 


which is originated from the one 
results therein and also in Ref. 
much larger than those to Cmm and C, 


mrj' 


As given in Eq. fl2U]) that rj is proportional to ^dCee, the upper bound on |(6 ,jCp| can 




59 


65| 


m 


be obtained by comparing with the current experimental sensitivities on 
Table HI where we have used Mh = 70GeV, Ma = 95GeV, Ms^ = 310 GeV, Ms^ = 90 GeV, 
Mjy = 200 GeV, and sg = 0.1. The strongest constraint is l^dCeel < 2.8 x 10“^, given by the 
target nucleus Xe. Finally, the contribution from Oj (proportional to Crfq) is much larger 
than that from {My)ee (proportional to Cmm) by a factor of 0(10“^), because the latter 
is greatly suppressed by the factor of me/v. On the other hand, lifting up \{My)ee\ to the 
average size of My ~ 10“^ eV will result in the excess of the 0z//9/3 decay rates that conflict 
with the observations. Table [I] also shows the maximum value of {My)ee for each nucleus. 


TABLE I: Constraints on I^^Ceel from for different nuclei as the targets. The corresponding 

limitation on \My\ee are also given. 


> Te,p(1025yr) 

Crin (lO^Syr)-! 

I'^dCeelmax 

|M,Ue(10-2)eV 

CERDA-1 (76 Qg) ^ 

2.1 

4.4 X 10"^ 

3.6 X 10-^ 

< 0.017 

KamLAND-Zen(736Xe) 

1.9 

8.3 X 10“*^ 

2.8 X 10“^ 

< 0.013 

NEMO-3(760Nd) [6^ 

0.0018 

2.9 X 10-7 

1.5 X 10-3 

< 0.072 

CUORICINO(730Te) [63j 

0.3 

2.3 X 10“® 

4.2 X 10-"^ 

< 0.02 

NEMO-3(®2Se) [64,^ 

0.036 

1.5 X 10“® 

1.5 X 10-3 

< 0.07 

NEMO-3(700Mo) 

0.11 

3.5 X 10-® 

5.6 X lO-'^ 

< 0.027 
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Finally, we end this section by mentioning that the role of the Z 2 symmetry is to make 
the dimension-7 operator O'j become the dominant contributions to in the LNV processes 
and Majorana neutrino masses. Note that as the quantum numbers of the doublet (y) 
and the singlet (s) scalars are the same as those in the Zee’s model 27|, the Majorana 
neutrino masses would be mainly generated by the corresponding one-loop diagrams related 
to the conventional Weinberg operator if the Z 2 symmetry is absent. However, with the Z 2 
symmetry, O'j is singled out at 1-loop level, while other LNV effective operators, especially 
the dimension-5 Weinberg operator, are much suppressed since they would be only induced 
by higher-loop diagrams. In this way, the Z 2 symmetry breaks the conventional effective 
operator ordering based on the scaling dimensions. Other LNV effects, like 0z//5/9 decays, 
would also change the leading modes accordingly. 


III. PHENOMENOLOGICAL CONSTRAINTS 


A. Electroweak Precision Tests 


As discussed previously, in order to have the two-loop neutrino masses in our model, 
the non-zero coupling constants A 5 and k are both required. The former splits the masses 
between 77° and A°, and the latter mixes the charged states and which carry different 
EW gauge quantum numbers. Both effects could change the values of the EW oblique S 


and T parameters. In particular, the T parameter shou 


d yield a stronger constraint on this 


model. The deviation of T from the SM is given by 13| 

I rr.2 „2 


AT = 






( 21 ) 


with the function F dehned by 


F = 


Ml + M2 


M^ 


m2 m2 

-" ^ log , 

X y y 


( 22 ) 


The value of becomes zero when -A- My, and it increases with the mass splitting 
among the new scalars. Note that AT has little to do with Di since there is neither mixing 
between Di and the SM leptons nor tree-level mass splitting among Di, while the deviation 
for the S parameter can also be ignored 6 ^. The formulae of Eq. (|2T|) is a general result 


10 













Mff=60GeV Mh=1QGqN 



FIG. 3: Plots of AT versus Ms^ with (a) Mh = 60GeV and (b) Mh = 70GeV, where blue, red, 
and green colors represent Ma — Mh = 25, 50, and 75 GeV, respectively, and solid (dashed) curve 
corresponds to sq = 0.1 (0), while the black dot is the benchmark point discussed in Sec. III.D. 


for the models with the mixings between the inert donblet and singlet scalars. The global 
htting resnlts constrain AT at 1.5 — 1.7a deviation by —0.02 < AT < 0.12 ^|. We show 
AT as a fnnction of Ms^ in Fig. [3l where we have nsed Ms^ = 310 GeV and sq = 0.1(0) 
along with (a) Mh = 60 and (b) 70 GeV. It is obvious that the numerical result of our 
model with se = 0.1 is approximately equal to a pure inert doublet model with se = 0. In 
general, AT goes up with increasing Ms 2 , and for a large value of Ma — Mh, the constraint 
on Ms 2 becomes stronger. The hgure also shows that Ma — Mh is limited to be less than 
75 GeV for Mh = 60 to 70 GeV and 90 < < 110 GeV. Finally, it should be noted that 

Ms 2 can not be too small, since there exists a lower bound on located within 70 to 
90 GeV from the LEP experiments. 


B. Dark Matter 


In this model, the lightest of the extra neutral particles: A^, and T )°2 3 could be a 

DM candidate, whose stability is guaranteed by the imposed Z 2 symmetry. In the following, 
we will focus on the case that DM is constituted solely by with a small charged scalar 
mixing S£, in which our DM would be very similar to that in the well-studied inert doublet 


model 


68 


70|. Furthermore, we concentrate on the low DM mass region with 50 GeV ^ 


Mh ^ 80 GeV 69|, [TOj, in which a large H^-A^ mass splitting can be allowed for the 
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generation of the right two-loop neutrino masses. In addition, the mass of should be 


higher than 90 GeV in order to escape the LEP bounds so that the co-annihilation 
channels, such as H^-A^ and would be strongly suppressed and thus ignored. 

We use the package micrOMEGAs 7l|] to accurately calculate the relic abundance in 


the above parameter space, including all possible annihilations and co-annihilations. When 
Mh approaches the half of the SM Higgs mass Mh/2 ~ 62.5 GeV 72, 73|, the Higgs res¬ 


onance in the s-channel would become prominent, which is characterized by the coupling 
\l controlling the trilinear vertex {XLv)hH^H^. However, in other regions, the DM anni¬ 
hilation cross section is dominated by the mode. Therefore, the correct DM relic 

abundance is achieved mainly by the balance of the and Higgs resonance channels. 


Fig. m shows the relevant parameter space in the Mh-\\A p 


served DM abundance 0.112 < Vtuh? < 0.128 at 3cr level 


ane, which can give the ob- 


49|. Note that when DM is 


heavier than 73 GeV, the WW^*'> channel would give a too large annihilation cross section 
to accommodate the DM relic abundance 69|, iTOj, which is omitted in Fig. IH 

The DM in this low mass region could be constrained by the DM direct detection 
experiments. Since we need a relatively large Higgs-mediation annihilation channel to gen¬ 


erate DM relic abundance, the Higgs exchange channel can also give rise to sizeab 


independent signals, with the corresponding DM-nucleon cross section as follows 69| : 


e spm- 




mt 


A/ 


^ -V 

dvr \MhMI) 


fm 


N ■ 


(23) 


Gurrently, the most stringent bound on the spin-independent DM-nucleon cross section is 
provided by the LUX experiment j^, with the minimum cross section of 7.6 x 10“^® cm^ for 
a DM mass of 33 GeV. It is shown in Fig. |4]that the LUX experiment has already probed 
some parameter space required by the DM relic abundance. Especially, the low DM mass 
region with Mh ^ 52 GeV is actually ruled out, as indicated by the shaded area in the plot. 
However, most parameter spaces are still allowed by LUX. 

C. Lepton Flavor Violation 


The current experimental constrain 
I —>■ /'y 3, 3) ~ ^ conversions 77 


s on LFV processes, such as the radiative decays 


80| . and three-lepton decays I hhh [Ml, ImI, are 


all dominated by one-loop diagrams with Z 2 odd particles inside. We take /i —)■ ey as an 
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FIG. 4: The parameter space for the correction DM relic abundance in the Mh — \)^l\ plane, with the 
blue band representing the region within 3a deviation of the cold dark matter relic abundance. The 
gray shaded area is excluded by the LUX experiment, and the black dot represents the benchmark 
point. This plot is calculated with Ma = 95 GeV and = 90 GeV. 


illustration because the current experimental upper bound < 5.7 x 10 

usually constrains a model in the most stringent way. In our model, we have 

r(/i"F ^ e^'j) 


Br(p^ —)■ e^y) = 


r(/i^ —>• e^uu) 

ScXp I ^ > 

/ . s/mC 


647rG|. 


le 


SqKsi + cIKs2 + 2^'^ 2^'^ 


(24) 


where the loop integrals and K'^ are dehned as 
2 z^ + 3z^ — 6z + 1 — 6z‘^ log z 


K. = 


6(1 - z)^M^ 


,K = 


z^ — 6z^ + 3z + 2 + 6z\ogz 

6(1 - z)*M^ 


(25) 


with 2 ; = Mfy/Ml- As is expected, the decay branching ratio is proportional to the coupling 

2 

constant combination J2iCifj,Cie ■ On the other hand, when the values of ^dCii' hxed 
by the neutrino masses, the only degree of freedom left is the size of ^d- When ^d is large, 
(ill should be suppressed, along with all the relevant LFV processes. This feature can 
be displayed in Fig. [5l where (w are expressed in the forms given in Eq. (fT8|) . with the 
unknown (eg satisfying the relation |.CdCee| ^ 10“"^, as well as se = 0.1, Mh = 70, Ma = 95, 
Ms^ = 310, Ms 2 = 90 and 120 and Md = 200 GeV. We hnd that the texture 7 c(d) yields 
the most stringent (weakest) constraint on (d, such that (d ^ 0.01 (0.002) is required for 
Ms 2 = 90 GeV. From another angle, if we set (d = 0.005 with Ta, we can predict that 


13 














M5^=90GeV 


M5^ = 120GeV 




FIG. 5; Br(// —)• 67 ) versus with (a) = 90GeV and (b) = 120 GeV, where the blue, 

red, green, and yellow curves correspond to the neutrino textures Ta, Tg, Tc, and Tg, respectively, 
while the black dot is the benchmark point. 

Br(/i — 67 ) = 10“^^ (3 X 10“^®) for = 90 (120) GeV, which might be measured by the 
next-generation experiments in the future. 


D. Numerical Results 

Based on the above constraints from the LFV processes, EW precision tests, direct 
searches of DM with the required relic abundance, we hnd a benchmark point from the 
allowed parameter space, given by: 


Mh = 70 GeV , Ma = 95 GeV , Mg, = 310 GeV , Mg, = 90 GeV , 

Md^ = Md^ = Md^ = Md = 200GeV , se = 0.1 , = 0.01, (26) 


and 


/ 




0.005 

0 . 

0 . 


0 . 

0.005 

0 . 


0 . 

0 . 

0.005 




, C 


^ 0.02 0.005 0.0022 ^ 

0.005 0.038 0.0061 
^0.0022 0.0061 0.0029^ 


(27) 


It is clear that the matrix ( corresponds to the neutrino texture Ta- We also plot this 
benchmark point by a black dot in Figs. [5^, [Sb, and 01 where the experimental results 
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from LFV processes, oblique parameters, and DM searches are well satisfied, respectively. 
Finally, from the benchmark point, one can obtain A 3 = 0.155, A 4 = —0.067, A 5 = —0.068, 
K, = 50.3 GeV, = 65.5 GeV, and Ms = 309 GeV, which are all small enough to ensure the 
self-consistence of the perturbation theory. 


IV. CONCLUSIONS 

We have tried to make the connection between neutrino physics and dark matter searches. 
In particular, we have emphasized that every effective operator, which violates the lepton 
number by two units, can give an equally good mechanism to generate Majorana neutrino 
masses. The problem lies in the fact that the new high-dimensional operators might be 
buried by the overwhelming effects from the conventional Weinberg operator which possess 
the smallest scaling dimension. One way to break this effective held theory ordering is to 
impose some symmetry which would protect the lightest neutral symmetry-protected states 
to be the DM particle. 

We have explicitly realized this connection by constructing a UV complete model with the 
Z 2 symmetry to generate the dimension-7 operator Oj = We have shown 

that the Majorana neutrino mass matrix structure and the leading 0z//9/9 decay contribution 
are closely related to O 7 . Especially, the neutrino mass matrix is predicted to be of the 
normal ordering due to the hierarchy in the charged lepton masses, and the 0z//5/9 decay 
rate can be large enough to be tested in the next-generation experiments. If we impose 
an additional GP symmetry in the lepton sector, we can even determine the form of the 
neutrino mass matrix completely. We have also focused on a specihc parameter region with 
a small mixing between charged scalars, and considered the constraints from the electroweak 
precision tests, dark matter searches, and LFV processes. 
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